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Effects of nonmagnetic impurity doping on an antiferromagnetic spin- 1/2 Heisenberg ladder 
system are studied by the quantum Monte Carlo method. A single nonmagnetic impurity induces 
a localized spin-1/2 moment accompanied by static and enhanced antiferromagnetic correlations 
around it. Small and finite concentration of impurities induces a remarkable change of mag- 
netic and thermodynamic properties with gapless excitations. It also shows rather sharp but 
continuous crossover around the concentration of about 4%. Above the crossover concentration, 
all the spins are strongly coupled participating in the enhanced and rather uniform power-law 
decay of the antiferromagnetic correlation. Below the crossover, each impurity forms an antifer- 
romagnetic cluster only weakly coupled each other. For random distribution of impurities, large 
Curie-like susceptibility accompanied with small residual entropy is obtained at low tempera- 
tures in agreement with recent experimental observation in Zn-doped SrCu203. Temperature 
dependence of antiferromagnetic susceptibility shows power- law- like but weaker divergence than 
the single chain antiferromagnet in the temperature range studied. 
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Spin-gapped Mott insulators have recently attracted 
much attention since the discovery of high-temperature 
superconductors with pseudo-spin-gap phase near the 
Mott insulator. A typical theoretical model for the spin- 
gapped Mott insulator is found in the antiferromagnetic 
Heisenberg (AFH) model with 2-leg ladder shape (two 
coupled spin chains) where a large spin gap is basic^ 
due to the dominant singlet formation on the rungs 
The spin gap is numerically estimated as A ~ 0.5Jl£ 
for the uniform exchange coupling J. Experimentally, 
the ladder model is believed to be relevant for systems 
such as (VO)2P207 and SrCu203. Recently, effects of 
Zn-doping on SrCu203 (Sr(Cui_xZnx)203 ) have been 
studied where Zn^+ ion substitutes Cx2+ {S = 1/2) and 
plays a role of nonmagnetic impurity.QQ-' Even small con- 
centration of nonmagnetic impurities causes a drastic 
change of magnetic properties with disappearance of the 
spin gap. At low temperatures below lOK, the Zn-doping 
induces anomaly of the susceptibility which is perhaps 
due to the antiferromagnetic (AF) transition. Above 
the anomaly of presumable Neel temperature, Zn-doping 
performed for the concentration 0.01 ^ x ^ 0.08 shows 
the Curie-like susceptibility which is consistent with the 
formation of free localized moments proportional to the 
number of impurities. On the other hand, the T-linear 
dependence of the specific heat is observed upon doping 
above Tn with the 7 value per spin similar to that of ID 
AFH chain (=2/3 J). The low temperature specific heat 
implies that the residual entropy just above Tn is incon- 
sistently smaller than that is expected from the free local- 
ized spins with the Curie law because J ^ lO^K or larger. 
Recently exact diagonalization of clusters and the varia- 
tional Monte Carlo calculations have been performed for 



the ladder model with regularly distributed nonmagnetic 
impurities .eP It suggests quick collapse of bulk gap struc- 
ture in the excitation spectra with sensitive enhancement 
of AF correlations for small impurity concentration con- 
sistently with experimental indications. This model has 
also been studied with the help of mapping to the ran- 
dom bond model and the appearance of the Curie-law 
has been suggested in three different regions of temper- 
atures.cl' 

In this paper, to understand these somewhat puzzling 
experimental results with sensitive dependence on the 
impurity concentration, depleted ladder system is stud- 
ied using an unbiased numerical method for large sys- 
tem sizes. We show results of Monte Carlo calcula- 
tions for the ladder model with nonmagnetic impuri- 
ties in both cases of regular and random distributions 
of impurities. In our results, large Curie-like suscepti- 
bility and small residual entropy are consistently under- 
stood within this model, where the overlap of AF clusters 
formed around impurities is substantial. We also show 
that rather sharp crossover exists around the impurity 
concentration Xc ~ 0.04. This is due to the fact that 
the overlap of AF clusters formed around the impurities 
becomes exponentially weak below this impurity concen- 
tration. To understand the low temperature properties 
for X < Xc, we point out the necessity to take account 
of three dimensional inter-ladder coupling to explain the 
experimental indications even above Tn. 

Our model is 2-leg ladder AFH model with periodic 
boundary condition given by the Hamiltonian, H = 
JJ2{ij) ^i' 1 where the summation is taken for all pairs 
of nearest neighbor sites. To take into account the non- 
magnetic impurity effect, we deplete the spins in the im- 
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purity sites. For the convenience of notation, we define 
an index of site as (x, y), where x{y) is the coordinate in 
the leg(rung)-direction of the ladder. For instance, for 
64x2-site system, x is from 1 to 64 and y is from 1 to 2. 
The site {x, y) belongs to A(B)-sublattice when x + y is 
even(odd) number. 

We adopt the quantum Monte Carlo (JDJvIC) method 
with an efficient loop cluster algorithm.lS'EP This algo- 
rithm together with the improved estimators reduces the 
statistical errors of the observables and the simulation 
time. We treat 64x2-site system with 0, 2, 4, 8, 10 
and 16 impurities (namely, impurity concentration x=0, 
0.0156, 0.0313, 0.0625, 0.0938, 0.125), and 60x2-site sys- 
tem with 10 impurities (a::=0.0833). The lowest temper- 
ature in our simulation is 1/120 in the scale of J. In 
terms of experiments in the Cu oxide ladders, this low- 
est temperature corresponds to the order of lOK which is 
indeed relevant when one wishes to compare with experi- 
mental results. Simulations are performed with 10"^-timc 
sweeps for equilibration and 5xl0^-time sweeps for mea- 
surements. The autocorrelation time of the loop cluster 
algorithm is typically a few sweeps. In case of the ran- 
dom distribution of impurities, we take typically average 
over 72 configurations of samples. 
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Fig. 1. l/T-dependence of the uniform susceptibility for 0(o), 
2{n), 4{0), 8(x), and 16(A) impurities on 64x2-site system. 
Lines A, B and C represent the fitting by the theoretical form. 
Other dashed lines are drawn as a guide to eye. The inset 
shows T-dependence of the uniform susceptibility for 64site AFH 
chain(») and ladder system with 16 impurities(A). 



First we investigate effects of nonmagnetic impuri- 
ties distributed regularly, where one impurity is lo- 
cated at site (1,1) and other impurities are put on two 
legs alternatingly with the same interval. In Fig.l, 
we show the temperature dependence of susceptibility 
X = T-^N~^Y^^ .{S^S^), where N, is the number of 
sites. This figure shows the qualitative difference be- 
tween 4- and 8-impurity. Without an impurity, because 
of the spin gap A, susceptibility becomes zero at low 
temperatures in accordance with the theoretical formEP 
X(T) ~ CoT-i/2exp(-A/r) (Line A), where Cq is a 
constant. With 2 and 4 impurities, within the tempera- 
ture range we studied, the susceptibility is well described 



"oyxiT) - CoT-i/2exp(-A/T) + Ci(n)/T (Lines B and 
C), where Ci{n) — n/ANg is the Curie constant which 
is proportional to the number of impurities n. The spin 
gap is estimated as A = 0.497 ±0.001 in the pure system 
while the parameter A is estimated as A = 0.526 ±0.001 
and A = 0.553 ± 0.001 for 2- and 4-impurity cases, re- 
spectively. At a; < Xc{^ 0.04), each impurity generates 
one spin- 1/2 moment in the sea of spin-gapped back- 
ground. Because the coupling between these spin mo- 
ments is weak, they behave as nearly free in this tempera- 
ture range. For example, for 2-impurity case at l/r=40, 
the two impurities are sufficiently well separated so that 
we can see effects of single-impurity doping around one of 
it. Each isolated impurity induces S'=l/2 spin. In addi- 
tion it generates static correlation whose saturated value 
at large r for sites ro=(l,2) and rj = {j,2) around the im- 
purity (1,1) shows an exponential dependence on |ro — fj| 
as (S'^(t=20)5;(0)) oc (-l)!'-""''^! •exp(-|ro-rj|/C) with 
^=2.91±0.15, which shows good agreement witLthe cor- 
relation length of pure ladder system, ^ ~ 3.19.lP There- 
fore it is concluded that an impurity induces 3—1/2 spin 
as well as static and enhanced AF correlation within the 
range of original correlation length of the gapped system. 
Similar condiisiDns have been reached by other theoret- 
ical studies .BEJ 

In contrast to this low impurity concentration, at 
X > Xc, the susceptibility does not show the Curie- 
law behavior in this regular distribution of impurities. 
As is shown in the inset of Fig.l, it is rather similar 
to that in the ID pure AFH with gapless spin excita- 
tion. (Note that the actual data suffer from finite-size 
effect at very low temperatures with vanishing suscep- 
tibility as T — > 0). Spins localized around impurities 
become strongly coupled above the crossover concentra- 
tion Xc ~ 0.04 in the temperature range studied. 
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Fig. 2. T-dependence of specific heat per spin for 0, 2, 4, 8, 10 
and 16 impurities. Line A shows the theoretical form for the 
pure ladder at low T. Line B represents the specific heat per site 
for ID AFH chain. Other lines are drawn as a guide to eye. 



In Fig. 2, temperature dependence of the specific 
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heat per spin is presented. For the pure ladder sys- 
tem, if we assume the momentum dependence of the 
low energy triplet dispersion as e(fc) = A-|-a|A: — 
TTp, with a being a constant, _low temperature (T <^ 
A) form of the specific heatEP is given by C(T) = 



correlation S[Q={n ,tt),t=Q) = N^^ ■ e 



'^/^. This 



form fits well to our QMC data for T <, 0.2, with the 
spin gap A estimated as 0.49 ± 0.04. At x < Xc{^ 0.04), 
temperature dependence of the specific heat looks simi- 
lar to the pure ladder case. For these dopings, as is seen 
from the data of the spin susceptibility, the systems are 
described by the free spins induced around the impuri- 
ties in our temperature range. Similarity of C to the un- 
doped system implies that the bulk gap structure in the 
excitation spectrum remains unchanged. At low enough 
temperatures, small interaction between these localized 
spins is expected to yield small peak in the specific heat 
to release the entropy, which is beyond the temperature 
range of our QMC simulations. 

As is seen in Fig. 2, at a; > Xc, the specific heat at 
low temperatures (T ^ 0.1) is enhanced and looks qual- 
itatively different from the case of low doping. The low 
temperature tails of the specific heat appears to follow 
T-linear dependence, as is pronouncedly visible for 16- 
impurity doping. The T-linear dependence is character- 
istic to the ID system with gapless excitations. In fact, 
for ID AFH, from the Bethe Ansatz calculation, low 
temperature specific heat is given as C = (2/3J)TE3^ 
per site. This value fits well for the T-linear part ob- 
served at a; > Xc. It implies that below a crossover 
temperature Tcr, all the spins become strongly coupled 
similarly to the ID AFH model. At x > Xc, Tcr in- 
creases more or less linearly with the doping concentra- 
tion, namely, Tcr ~ 0.04, 0.06 and 0.09 for n = 8, 10 and 
16, respectively. In fact, when the entropy of n localized 
spins starts released below Tcr with T-linear specific heat 
C - (2/3J)T, Tcr is estimated to be Tcr ~ (3Ja;/2)ln2. 
This is consistent with the numerical data for x > Xc- Al- 
though the exponent of the specific heat a defined from 
C oc T" at asymptotically low temperatures may be dif- 
ferent from 1 for the randomly distributed impurities, the 
residual entropy in our temperature range may be similar 
between regular and random distribution of impurities. 
Therefore, the specific heat at x > Xc is consistent with 
the experimental observation in Zn-doped SrCu203. The 
peak in the specific heat observed below Tcr is consistent 
with temperature dependence of where localized spins 
become strongly coupled below Tcr. From the calculated 
specific heat, for Zn doped ladder system, we predict 
that the specific heat linear in T should decrease at Tcr 
when we increase temperature further before the expo- 
nential growth due to the bulk spin gap appears. In the 
experiments, the linear specific heat is observed even at 
the concentration as low as a: 0.01 around 10 K. At 
least in the temperature of the order of lOK, we could 
not identify the linear specific heat for x < Xc- This dis- 
crepancy may be due to the three dimensional coupling. 
In fact the AF order appears to be set in around the 
similar temperature. 

In Fig. 3, we plot T-dependence of equal time AF spin 



and AF spin susceptibility x(Q=(7i',7r),Ci;=:0) = N^^ j^^^ dr 
X;,.j e-*'3-('''-'^^)(S'f(r)S';(0)) for both regular and ran- 
dom impurity distributions. 

For the regularly distributed case, in contrast to the 
4-impurity case, 8-impurity case shows good coincidence 
with ID AFH for T ^0.1, while 5(Q=(7r,7r),r=0) shows 
similar T-dependence as 0-impurity case for 4-impurity 
case. It turns out that a,t x < Xc, spin gap nature of the 
ladder system is not changed, and the correlation decays 
exponentially in the temperature range studied. The im- 
purity spins become coupled with presumable power-law 
decay of correlations only below Tcr. At a; < Xc, Tcr de- 
creases exponentially fast with decreasing x and it disap- 
pears from experimentally accessible temperature range 
quickly because the effective exchange coupling between 
neighboring localized spins confined around the impuri- 
ties should decrease exponentially with the decrease in 
the doping concentration x. However at a; > a;c, spin 
correlations follow the power law decay, which again in- 
dicates that all the spins become strongly coupled as in 
the ID AFH chain. In fact in this regular distribution of 
impurities, spin correlations in real space have quantita- 
tively similar values to the case of ID AFH chain with 
the same power oc 1/r. An important point is that in 
this high-doping region, not only the localized spins con- 
fined around impurities but all the spins in the system 
take part in this enhancement of spin correlation. This is 
consistent with the results of the magnetic susceptibility 
and the specific heat. 
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Fig. 3. Equal time spin correlation 5(Q=(7r,7r),r=0) and an- 
tiferromagnetic suscpetibility x(Q=(''''i''")i'^=0) a-re plotted for 
0(x), 4(o), 8(A) impurities with regular distributions while bold 
curves show the case of 128-site pure AFH chain. Filled sysm- 
bols show the case of 4 (circle) and 8 (triangle) impurities with 
random distributions. 



Next we discuss the case of random impurity distribu- 
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Fig. 4. T-dependence of inverse uniform susceptibility for 2(n), 
4(0), 8(x), 12(») and 16(A) impurities with random distribu- 
tion on 64x2-site system. The inset shows the Curie constants 
obtained from our calculation(o) and the experiment(»). 



tion. Residual entropy is expected to be similar to the 
case of regular impurity distribution in our temperature 
range. However, as we see in Fig. 4, the uniform sus- 
ceptibility has qualitatively different behavior. For all 
impurity concentrations, even at x > Xc, the Curie- like 
susceptibility is observed, In fact, from the arguments 
by Sigrist and Furusakijj) at T ^ A (spin gap of pure 
ladder system), Curie behavior of uniform susceptibility 
is expected with two different Curie constants a;/4 and 
x/12 with a crossover between these two Curie behaviors 
around T=T*. Here T* is the order of average coupling 
between impurity spins. Below T* all impurity spins 
couple with each other and the system shows the spin 
moment due to the fluctuating imbalance of the number 
of impurities on A and B sublattices. In the relevant tem- 
perature range in the experiment ~ 0.0U(< A), our cal- 
culated Curie constants seem to lie between x/4 and s/12 
and is consistent with the experimental observation. ffl In 
terms of the argument by Sigrist and Furusakijj^ it im- 
plies that there is a wide intermediate temperature re- 
gion between these two Curie behaviors. 

For the case of randomly distributed impuri- 
ties in Fig. 3, the growths of S{Q={TT,Tr),T=0) and 
x(Q=(7r,7r),u;=0) at low temperatures are suppressed 
both from the cases of regularly distributed case and 
the ID AFH case. S'(Q=(7r,7r),T=0) shows very weak 
temperature dependence in the temperature range stud- 
ied. x(Q=(7r,7r),a'=0) appears to follow approximately 
power-like divergence cx T^" (a ^ 0.7) with weaker diver- 
gence than the case of ID AFH (cx T^^). Although the 
AF correlation is clearly enhanced from the pure ladder, 
its enhancement is not as large as the ID AFH case. This 
numerical result does not exclude stronger enhancement 
of AF correlation at lower temperatures. However in the 
experimentally relevant temperature range T ~ O.OU, 
this-,result has shown that the argument of Nagaosa et 
sl\S3 is not applicable because they predict a more sin- 
gular divergence of x((5=(7r,7r),w=0) than ID AFH. 

In summary, effects of nonmagnetic impurities on a 
spin ladder are investigated by the quantum Monte Carlo 



method with loop cluster algorithm. Single impurity 
causes the formation of a spin- 1/2 localized spin con- 
fined around the impurity site. This impurity spin gen- 
erates a cluster of the size ^ ~ 3 around it with static 
AF correlation, which is the origin of the enhancement 
of the AF correlation. From the results from regularly 
and randomly distributed impurities, a rather sharp but 
continuous crossover concentration (xc ~ 0.04) is identi- 
fied. Above Xc, all the spins are strongly coupled with 
rather uniform and power-law decay of the AF corre- 
lation. Below Xc, the impurity spins are only weakly 
coupled and is expected to give rise to power-law de- 
cay of the correlations between impurity spins only at 
sufficiently low temperatures. Above Xc our results are 
consistent with experimental results in terms of two as- 
pects. One is a large Curie-like susceptibility with the 
Curie constant x. The other is a small residual en- 
tropy and T-linear specific heat with 7 similar to the 
per spin value of ID AFH chain. However, it seems to 
be hard to keep the same behavior in the temperature 
range of O.OIJ at x = 0.01. We speculate that the three 
dimensional coupling plays an important role in keeping 
the similar behavior even at x~ 0.01. The AF correla- 
tion is strongly enhanced as compared to the pure ladder. 
However, its enhancement is not as large as the ID AFH 
chain. 
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